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We investigate a trajectory for the Wilson flow in the theory space. For this purpose, we deter¬ 
mine the coefficient of the plaquette and rectangular terms in the action for the configurations 
defined by the solution of the Wilson flow. The demon method regarded as one of the inverse 
Monte Carlo methods is used for the determination of them. Starting from the conventional 
Wilson plaquette action of quenched QCD, we find that the coefficient of the plaquette grows 
while that of the rectangular tends to negative with the development of the flow as the known 
improved actions. We also find that the trajectory forms a straight line in the two-coupling 
theory space. 
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1 Introduction 


The lattice gauge theory is one of the most powerful tool for investigating the quantum 
field theory. However, it is not suitable for studying quantities which are sensitive to the 
cutoff: the topological charge, energy momentum tensor and so on. The Wilson flow is con¬ 


sidered to be a remedy for such issues [1R9] • It is regarded as a continuous stout smearing and 
one can measure such observables even on the discretized spacetime [3,llH- In Ref. lfl|. M. 
Liischer has studied the effective action for the configuration at finite t = t/a 2 in the context 
of the transformation in the field space and given a analytic form of it. The purpose of this 
paper is to determine it numerically by the so-called demon method. 

The demon method has been proposed by M. Creutz [uj and improved by M. Hasenbusch 
et al. |13?], which enables us to determine effective couplings from a given configuration. 
There has been various applications of this method. For example, investigating the couplings 
induced by the fermionic determinant (mJ, studying the renormalization group of the SU(3 ) 
lattice gauge theory flil. a nd describing the SU (3) lattice gauge theory in terms of effective 
Polyakov-loop models 16] . Using the method, one can investigate effective actions without 
any perturbative approximation. 

A similar work to ours has been done by QCD-TARO collaboration [ 17 I 18 ]. They have 
studied the renormalization trajectory for quenched QCD in the two-coupling theory space 
performing the blocking of the link variables [h]] and using the Schwinger-Dyson equation 
method [20J to compute the effective coupling of the action for the blocked configurations. As 
mentioned above, since the Wilson flow is regarded as a continuous smearing, we compare 
our effective action with the ones investigated in Ref. [wl, 18] . As we will see later, our 
effective action shows the same tendency as the known improved action, however it has 
different pictures. 


In this work, we determine a trajectory for the Wilson flow in the theory space using 
the demon method. We choose /3 = 6.0 Wilson plaquette action for configuration generation 
and apply the Wilson flow. The range of the flow time is taken to 0 < V&t < 1.3. In order 
to evaluate finite size effects, our lattice sizes are taken to T 4 = (T/a) 4 = 4 4 , 8 4 and 16 4 . 
We find that the coefficient of the plaquette grows while that of the rectangular tends to 
negative with the flow time as the known improved actions. We also find that the trajectory 
forms a straight line in the two-coupling theory space. 

This paper is organized as follows. In Sec. 0J we briefly review the demon method and 
explain our strategy to obtain the effective action. A numerical setup and results are shown 
in Sec. [3l Finally, a conclusion, discussions and future perspectives are given in Sec. 01 
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2 The demon method and our strategy 


In this section, we review the demon method 12, |13j and how we utilize it for the 
determination of the effective action for configurations which are continuously smeared by 
the Wilson flow equation l|: 


V M (x;t) = -g%{d XiM S w (V)}V M (x;t), V M (x;t) | t =o = U^x), 


( 1 ) 


where S\y(17) is defined by 


1 _ ReTr W™(x) 


Swto =p y 

x,n<u L 

W^(x) = Uf,(x)U u (x + n)Ul(x + v)Ul{x). 


B= A 
$ 


( 2 ) 


For the definition of the effective action, let us consider an expectation value of any 
operator 0[Vt(U)\ 

(0)t = \ J DU 0[Vt(U )] e-^wTl, (3) 

where 

Z = I DU e ~ s ^ u \ (4) 


and Vt = {V^(x\ t)}. The effective action is defined by inserting the delta function and 
changing the integration variables from U to Vf as follows. 


where 


(0)t 




DV t S(Vt - V t {U )) 


DUDV t 8 (V l - V t (U)) 0[V t ] 
[ DV t 0[V t \ e~ Ses[Vt \ 


0[V t (U)\ e“ 5w ^ 
e -S w [u} 


(5) 


Zett = 


I 


DV t 


( 6 ) 


S'gfj can be regarded as the “improved” action since the discretization effects for observables 
are reduced at finite t. 
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2.1 The demon method 

Here, we briefly review the demon method which enables us to determine the couplings 
of the action from a given configuration. 

Let us consider a configuration whose distribution is the Boltzmann weight 

P(U) oc e~W u \ (7) 

where /3 is a priori unknown coupling and to be determined by the demon method. 

Now we introduce an extra degree of freedom “demon” to the system and consider a 
microcanonical partition function of a joint system, 

■^mic = EE S(S[U]+E d -Eo), (8) 

U E d 

where Eq is an initially determined total energy of the joint system, and E d is the energy 
carried by the demon. The demon energy is restricted within the range 77 m in < E d < E max 
which can be chosen suitably. Hereafter we set E m \ n = —E max for simplicity. Starting from a 
given configuration, we update the system keeping 5[f7] + Ed constant. The updating process 
is described as follows. First, a new configuration U' is proposed. To keep the total energy 
constant, the new demon energy is given by 

E’ d = E d - spy'] + spy], (9) 

When E' d is in the allowed region [—E majt _:E max ], the new configuration and the new demon 
energy are accepted, otherwise they are rejected and the configuration and the demon 
energy remain to be same. If the probability of the change is symmetric in {U, E and 
{U',E' d }, a generated sequence of combined configurations is distributed according to a 
uniform distribution. 

When the degrees of freedom of U are sufficiently large, a given configuration behaves as 
a heat bath to thermalize the demon. Therefore standard statistical mechanics arguments 
show that Ed is distributed with 

P(E d ) cx e~P Ed , (10) 

in the thermodynamic limit. The average of E d is given by 

1 rEmax -| 

(E d ) = - dE d E d e~P Ed = - 

J -Emax ' 

where 

/ -Emax 

dE d (12) 

' E/max 

Finally, /? is obtained by measuring (E d ) and solving Eq. (fTTTh 


1 - 


(3E n 


tanh(/3F/ max ) 


( 11 ) 
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The extension to a multi-coupling system is straightforward. Suppose that the action is 
parametrized by 


i 


(13) 


where Si[U] is the interaction term and Pi is the corresponding coupling. Now we introduce 
demons E\ for each coupling. During the microcanonical update, a proposed configuration 
and demon energies are accepted only if all demon energies are in the allowed region. Here we 
take the common allowed region [—i? max : E nmx \ to all couplings for simplicity, Pi is obtained 
by solving 

= S 


1 - 


PiEv 


tanh(/3j£' r) 


(14) 


2.2 Our strategy 

In this subsection, we explain our strategy to obtain the effective action defined by Eq. (jHJ). 
We assume that S e g can be parametrized by the form of Eq. m, where Si[U] denotes several 
types of Wilson loops such as the plaquette, rectangular, chair, sofa, etc. The simulation is 
implemented as the following steps. 

(1) Evolve a given configuration by integrating the Wilson flow equation Eq. (JTj) from 
0 to t. 

(2) Perform the microcanonical updates for the joint system using the configuration 
at t according to the procedure described in Subsec. 12.11 

(3) After a sufficiently large number of microcanonical updates, take an average of the 
demon energy and obtain /3j by solving Eq. m- 

(4) Replace the configuration at t = 0 by a new statistically independent one. This is 
done by the HMC updates of the original system at t — 0. In this step the demon 
degrees of freedom are frozen. 

(5) Go back to the step (1) with the new configuration. The initial demons’ energies 
in the step (2) are given by the averages of the previous run. 

The flowchart of the above algorithm is shown schematically in Fig. [D By taking averages of 
PiS obtained for each gauge configuration, we finally determine the couplings of the effective 
action. The step (4) and (5) are needed to suppress the systematic errors coming from the 
finiteness of the volume 11311 . 
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(4) HMC HMC 

c TT , u /^plaq = 6.0 dplaq = 6.0 

{U^x)} — - *{Up{x)} —►■■■ 


1 Demon method 


(1) Wilson flow with 
i r flow time t 

' (2) (3) Demon method 


■ Demon method 


Wilson flow with 
flow time t 


' Demon method 


dplaq(0) 



/Wo) 

Averaging 

dplaq(0) 

drcct(0) 



drect (0) 


drect (0) 


dplaq (t) 


dplaq (t) 

Averaging 

dplaq (0 

drect {t) 


drect (c 


drect {t) 


Fig. 1 Schematic picture of our strategy. Bold texts with the number correspond to the 
explanation in the text. 


3 Numerical setup and results 


3.1 Numerical setup 


Here, we introduce our numerical setup. We employ the Wilson plaquette action defined 
by Eq. (J2]) with /? = 6.0 for configuration generations. Configurations are generated by the 
HMC algorithm. Our lattice sizes are (L/a) 4 =L 4 = 4 4 , 8 4 , and 16 4 . 

In this work, we implement the demon method with the truncation 


S eS (U) = dplaq ^ 



iReTr W ^ 1 



+ drect ^ ^ 



^ReTr R 


1x2 

fiv 



where R l ^ 2 (x) is the rectangular loop defined by 


= U fi (x)U t[ (x + + 2[i)u] i (x + v + fx)uj.(x + v)Ul(x). 


(15) 


(16) 


Each gauge conhguration is separated by 10 trajectories of the HMC. The total number of 
measurements is 500, 500 and 250 for L 4 = 4 4 , 8 4 and 16 4 respectively. The autocorrelation 
times defined by the average of the plaquette of each configuration are around 4.6, 11 and 
11 for I 4 = 4 4 , 8 4 and 16 4 . 

Our measurements are performed as follows. All of the initial demon energies are set to 
be 0 and E max is set to be 5.0. We perform 10000 microcanonical updates of the joint system 
with a given trial configuration for the thermalization of the demon energies. First 2000 
updates are discarded and the averages of the demon energies over the last 8000 updates 
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are taken. These averages are used as the initial demon energies for the next microcanonical 
updates. Following steps are repeated during the measurements. After a replacement of the 
gauge configuration, we always perform the 10000, 5000 and 5000 microcanonical updates 
for L 4 = 4 4 , 8 4 and 16 4 . First 2000, 1000 and 1000 updates are discarded and we take the 
averages of the demon energies over the last 8000, 4000 and 4000 updates for L 4 = 4 4 , 8 4 
and 16 4 respectively. The above procedures are implemented for each flow time. The range 
of the flow time is taken to 0 < Vst < 1.3. 


3.2 Results 


We show the results of the flow time dependence of the effective action in this subsection. 
The values of /3 p i aq and /3 rec t for 0 < t < 0.20 are listed in Table. CD We examine the consis¬ 
tency between the effective action for t = 0 and the initial action. Since our initial couplings 
(/^plaq; /^rect) = (6.0, 0.0) are reconstructed within 2 sigma, the demon method works well. 

The results show that the value of (3 p i aq grows with t (Fig. [2D- We expect that this is 
because the Wilson flow has a smoothing effect on the gauge field and lowers the value of 
the Wilson action as t increases 10]. 

On the other hand, the value of /3 rec t tends to a negative region and decreases with t (Fig. 
[2D. This fact seems reasonable since the negativeness of is a common feature of the known 


improved actions such as the Symanzik 21-23], Iwasaki 24, [25] and DBW2 action [18|. Note 


that our scheme is not based on any perturbative analysis and therefore defines the non- 
perturbatively “improved” action. Moreover, the action may be systematically improved by 
adding conceivable Wilson loops to the ansatz of S e g. 

/3piaq and /3 rec t are plotted in Fig. [21 where the horizontal axis indicates Vst, which is the 
effective range of smearing the link variables. Since the values of /3 p i aq and /3 rec t for each L 
at a fixed t coincide with each other, finite volume effects turn out to be irrelevant. 

/3piaq and /3 rec t for L 4 = 16 4 are fitted in the forms 


W«) = ^ 8B ‘. Aect(t) = C( 1 - 


8 Dt\ 


(17) 


yielding the numerical results 


A = 6.040(47), B = 1.906(7), C = 1.010(47), D = 1.911(35), (18) 

and plotted with the curved lines respectively in Fig. [21 This result shows that the t 
dependence of the exponent is common to these two couplings within the errors. 

We also plot the flow of the effective action in the two-coupling theory space in Fig. [3l 
Fig. Q] shows the plot around the origin of Fig. [3] Reflecting the fact that two couplings have 
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L 4 = 4 4 

L 4 

= 8 4 

L 4 = 

16 4 

t 

/3plaq 

P rect 

/^plaq 

A’ect 

ft plaq 

A’ect 

0.00 

5.972(23) 

0.009(9) 

6.072(38) 

-0.024(13) 

5.932(66) 

0.030(24) 

0.02 

8.116(44) 

-0.363(12) 

8.368(75) 

-0.387(17) 

8.087(58) 

-0.329(25) 

0.04 

10.856(59) 

-0.837(17) 

11.334(49) 

-0.895(20) 

11.088(70) 

-0.817(30) 

0.06 

14.741(83) 

-1.490(19) 

15.512(63) 

-1.602(24) 

15.120(141) 

-1.440(69) 

0.08 

19.789(86) 

-2.329(21) 

21.204(82) 

-2.565(30) 

20.749(148) 

-2.408(41) 

0.10 

26.671(189) -3.477(39) 

28.637(200) 

-3.791(72) 

28.660(243) 

-3.705(56) 

0.12 

35.881(190) -5.002(34) 

37.32(117) 

-4.622(701) 

37.622(841) 

-5.045(347) 

0.14 

46.793(619) -6.682(152) 

51.208(848) 

-7.207(523) 

51.535(289) 

-7.692(75) 

0.16 

64.23(600) -9.635(127) 

66.03(287) 

-8.48(205) 

68.986(557) 

-10.711(132) 

0.18 

81.79(260) -11.88(101) 

90.41(251) 

-12.86(185) 

90.930(801) 

-14.408(207) 

0.20 

112.06(373) -16.83(101) 

122.56(99) 

-19.663(216) 

106.69(776) 

-11.92(613) 


Table 1 The values of (3 p i aq and /3 rec t at 0 < t < 0.20 for L 4 = 4 4 , 8 4 and 16 4 . 


the common exponent as the function of t, the flow of the effective action forms a non-trivial 
straight line (Fig. EJ): 

/dplaq + 5.931 (31)/?rect = 6.216(121), (19) 

determined by a numerical fitting for L 4 = 16 4 with y 2 /d.o.f. = 0.857(371). 

4 Conclusion and Discussions 

We have investigated a trajectory for the Wilson flow in the theory space by the demon 
method. The effective action is truncated such that it has the plaquette and rectangular 
terms. We have measured the flow time dependence of d p i aq and /3 rec t, and found that dplaq 
increases while /3 rec t tends to the negative region with t. As shown in Fig. [3] and Fig. [4l we 
have found that the trajectory forms the straight line in two-coupling theory space, although 
we do not know whether this is just an accident or not. Since the value of /3 has been fixed 
throughout this work, it should be examined whether this fact is universal to all values of /3. 
A trajectory in a wider parameter space has to be investigated by adding more Wilson loops 
to the ansatz as a next step of this work. It is also interesting to study how the trajectory is 
changed when we employ other actions for the gradient flow. 

As noted in Sec. [TJ one can measure the quantities which are sensitive to the cutoff by 
using the Wilson flow. From this fact, our effective action is expected to be an “improved” 




















v'sl 


Fig. 2 The values of /3 p i a q and (3 rec t as a function of \/8£. Two dashed horizontal lines are 
at the initial value of /3 p i aq = 6.0 and /3 rec t = 0 respectively. Two curved lines are the fit lines 
for I/ 4 = 16 4 . The details of the fit are described in the text. 


action which reduces the discretization effects. Actually, we have found that the flow time 
dependence of dpi an and drect shows the same tendency as in the already known improved 


actions 


18 


2ll425||. In order to confirm this, the cutoff sensitive quantities, such as the 


topological charge and the energy momentum tensor, should be studied using our effective 
action. It is an interesting question whether such quantities are well-measured by tuning 
only two parameters in the action. The restoration of the rotational symmetry should also 
be checke d by measuring, for example, the difference between the on-axis and off-axis Wilson 
loops (171, lla). The advantages of our effective action are as follows. First, S e ff itself is totally 
well-defined without any ambiguity such as truncations in blocking and projections of the 
link variables. Second, it does not require any perturbative analysis. Finally, we have one 
tunable parameter: the flow time. 

Related to the improvement of the lattice action, lattice artifacts could be reduced by 
the Wilson flow due to its smearing property. There is a well-known spurious UV fixed 
point caused by lattice artifacts when we have the adjoint plaquette term or the fermionic 
determinant 


14j . The demon method combined with the Wilson flow, which is explained 


m 


Subsec. 12.21 can also be applied to this issue. Namely, one can check the issue by measuring the 
coupling of the adjoint plaquette term of configurations generated by the fundamental-adjoint 
mixed action. 
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ft plaq 


Fig. 3 The flow of the effective action in two-coupling theory space. The value of t 


becomes larger along the arrow. The doted line represen 
Other lines represent the Symanzik 


21 


23| , Iwasaki 


24 


s our numerical fit for L 4 = 16 4 


25 1 and DBW2 action 18|. We 


take — A-ect//3piaq = 0.05, 0.09073 and 0.11481 for the Symanzik, Iwasaki and DBW2 action 
respectively. 


A 


rect 



/^plaq 


Fig. 4 The same figure as Fig. [3] but around the origin. 
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Since the Wilson flow is similar to a continuous “block spin transformation”, it may 
be used to define the exact renormalization group of the lattice gauge theories. However, 
the trajectory for the Wilson flow in the theory space travels in the opposite direction to 
the renormalization trajectory investigated by QCD-TARO collaboration 17, |l8|. Therefore, 


our effective action cannot be regarded as an action obtained after a renormalization group 
transformation. This is because the Wilson flow itself is just a “blocking” and does not 
contain a rescaling, which is necessary for the renormalization group transformation. Just 
replacing the blocking procedure 19] used in Ref. [l3, 18] by the Wilson flow is considered 
as a way to construct a renormalization group scheme via the Wilson flow. This is also an 
interesting subject related to our work. 

Finally, we comment on two types of systematic errors in this work. 

One stems from the truncation of the effective action. This can be reduced by adding 
possible Wilson loops to our truncated action. Since the configurations at large t picks up 
the information of the original configurations at widely separated points, larger Wilson loops 
should be required with increasing the flow time. Therefore our truncation may be reasonable 
only in the small t regime. 

Another comes from the demon method itself. As noted in Subsec. 13.21 the Wilson flow 
lowers the value of the Wilson action as t increases. Since the gauge configurations are 
changed randomly during the microcanonical updates, most of attempts raise the value of 
Sw an d are rejected. Therefore it is hard to determine a large-valued (3i precisely in the 
demon method. This error is unavoidable as long as using the demon method, however there 
is another way to determine the couplings from a given configuration using the Schwinger- 
Dyson equation method [20 ]. 


Acknowledgements 

We thank Hidenori Fukaya, Hideo Matsufuru and Tetsuya Onogi for fruitful discussions and useful comments. 
We also thank Hiroshi Suzuki for his warm encouragements. Guido Cossu, Hideo Matsufuru, Jun-ichi Noaki 
and Satoru Ueda gave us gen erous support for the setup of the numerical calculations. This work is in part 
based on Bridge++ code [26} (http://bridge.kek.jp/Lattice-code/) and Iroiro+-1- code {27}. The numerical 
calculations were carried out on SR16000 at YITP in Kyoto University. R.Y. is supported in part by the 
Grand-in-Aid of the Japanese Ministry of Education (No. 15J01081). 

References 

[1] Martin Luscher, JHEP, 08, 071, [Erratum: JHEP03,092(2014)] (2010), arXiv:1006.4518. 

[2] Martin Luscher, PoS, LATTICE2010, 015 (2010), arXiv: 1009.5877. 

[3] Martin Luscher and Peter Weisz, JHEP, 02, 051 (2011), arXiv:1101.0963. 

[4] Szabolcs Borsanyi et al., JHEP, 09, 010 (2012), arXiv:1203.4469. 

[5] Szabolcs Borsanyi, Stephan Durr, Zoltan Fodor, Sandor D. Katz, Stefan Krieg, Thorsten Kurth, Simon Mages, 
Andreas Schafer, and Kalman K. Szabo (2012), arXiv:1205.0781. 

[6] Zoltan Fodor, Kieran Holland, Julius Kuti, Daniel Nogradi, and Chik Him Wong, JHEP, 11, 007 (2012), 
arXiv:1208.1051. 


11 









[7] Zoltan Fodor, Kieran Holland, Julius Kuti, Daniel Nogradi, and Chik Him Wong, PoS, LATTICE2012, 050 
(2012), arXiv:1211.3247. 

[8] Patrick Fritzsch and Alberto Ramos, JHEP, 1310, 008 (2013), arXiv:1301.4388. 

[9] Martin Luscher, JHEP, 1304, 123 (2013), arXiv: 1302.5246. 

[10] Martin Luscher, Commun.Math.Phys., 293, 899-919 (2010), arXiv:0907.5491. 

[11] Masayuki Asakawa, Tetsuo Hatsuda, Etsuko Itou, Masakiyo Kitazawa, and Hiroshi Suzuki, Phys. Rev., D90(l), 
011501 (2014), arXiv: 1312.7492. 

[12] Michael Creutz, Phys.Rev.Lett., 50, 1411 (1983). 

[13] M. Hasenbusch, K. Pinn, and C. Wieczerkowski, Phys.Lett., B338, 308-312 (1994), arXiv:hep-lat/9406019. 

[14] T. Blum, Carleton E. DeTar, Urs M. Heller, Leo Karkkainen, K. Rummukainen, and D. Toussaint, Nucl. Phys., 
B442, 301-316 (1995), arXiv:hep-lat/9412038. 

[15] T. Takaishi, Mod. Phys. Lett., A10, 503-514 (1995). 

[16] Christian Wozar, Tobias Kaestner, Andreas Wipf, and Thomas Heinzl, Phys. Rev., D76, 085004 (2007), 
arXiv:0704.2570. 

[17] P. de Forcrand et al. (1998), arXiv:hep-lat/9806008. 

[18] P. de Forcrand, M. Garcia Perez, T. Hashimoto, S. Hioki, H. Matsufuru, O. Miyamura, A. Nakamura, I. O. 
Stamatescu, T. Takaishi, and T. Umeda, Nucl. Phys., B577, 263-278 (2000), arXiv:hep-lat/9911033. 

[19] R. H. Swendsen, Phys. Rev. Lett., 42, 859-861 (1979). 

[20] Antonio Gonzalez-Arroyo and J. Salas, Phys. Lett., B214, 418 (1988). 

[21] P. Weisz, Nucl. Phys., B212, 1 (1983). 

[22] P. Weisz and R. Wohlert, Nucl. Phys., B236, 397, [Erratum: Nucl. Phys.B247,544(1984)] (1984). 

[23] M. Luscher and P. Weisz, Commun. Math. Phys., 97, 59, [Erratum: Commun. Math. Phys.98,433(1985)] (1985). 

[24] Y. Iwasaki (1983). 

[25] Y. Iwasaki, Nucl. Phys., B258, 141 156 (1985). 

[26] S. Ueda, S. Aoki, T. Aoyama, K. Kanaya, H. Matsufuru, S. Motoki, Y. Namekawa, H. Nemura, Y. Taniguchi, 
and N. Ukita, PoS, LATTICE2013, 412 (2014). 

[27] Guido Cossu, Jun Noaki, Shoji Hashimoto, Takashi Kaneko, Hidenori Fukaya, Peter A. Boyle, and Jun Doi, 
JLQCD IroIro++ lattice code on BG/Q, In Proceedings, 31st International Symposium, on Lattice Field Theory 
(Lattice 2013) (2013), arXiv:1311.0084. 


12 


